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Abstract 

• ^ Wc study the abelian sandpile model on a random binary tree. 

^ Using a transfer matrix approach introduced by Dhar & Majumdar, 

^ we prove exponential decay of correlations, and in a small supercrit- 

ical region (i.e., where the branching process survives with positive 
probability) exponential decay of avalanche sizes. This shows a phase 
transition phenomenon between exponential decay and power law de- 
cay of avalanche sizes. Our main technical tools are: (1) A recursion 
for the ratio between the numbers of weakly and strongly allowed con- 
figurations which is proved to have a well-defined stochastic solution; 
(2) quenched and annealed estimates of the eigenvalues of a product 
of n random transfer matrices. 
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1 Introduction 



The abelian sandpile model (ASM) is a thoroughly studied model both in 
the physics and in the mathematics literature see e.g. [HI [121 El [15] for 
recent review papers on the subject. In physics, it serves as a paradig- 
matic model of self-organized criticality (SOC). SOC is usually referred to 
as the phenomenon that the model exhibits power law decay of correlations 
or avalanche sizes, without fine-tuning any external parameters such as tem- 
perature or magnetic field. In mathematics, the ASM is connected to several 
combinatorial objects such as spanning trees, graph-orientations, dimers, and 
it has an interesting abelian group structure. 

The ASM has been studied on the Bethe lattice (i.e., the rootless binary 
tree) in |5] . Via a recursive analysis, based on a transfer matrix method, the 
authors in [5] arrive at exact expressions of various quantities of interest, such 
as the single height distribution, correlation functions of height variables, and 
avalanche size distribution. 

There are various motivations to consider the ASM on random graphs. As 
an example, we mention integrate-and-fire models in neuroscience (see e.g. 
[in] and related papers), where the connections between neurons are updated 
after a neuron has fired. The typical connection structure of a network of 
firing neurons is therefore generically not translation invariant, and time 
dependent. As a first approximation, one can quench the randomness of 
the connection graph and study the firing of neurons on the derived random 
graph. So far, the ASM model has been studied on small world graphs from 
the physicist's perspective, using a renormalization group approach [8j. In 
the mathematics literature, there are recent studies on so-called "cactus" 
graphs [13]. 

In this paper, we start this study of the ASM on random graphs with the 
ASM on a random tree, for the sake of simplicity chosen to be a realization 
of a binary branching process with branching probability p e [0, 1]. We use 
the transfer-matrix method of [5] to express relevant quantities such as the 
correlation of height variables and the avalanche size distribution in terms 
of the eigenvalues of an ad hoc product of random matrices. This is the 
fundamental difference between the Bethe lattice case and the random tree, 
namely the fact that the transfer matrices depend randomly on the vertices 
and instead of having to deal with the n-th power of a simple two by two 
matrix, one has to control the product of n random matrices. 

The crucial quantity entering the transfer matrices is the so-called charac- 
teristic ratio, which is the ratio between the numbers of weakly and strongly 
allowed configurations. This ratio is equal to 1 in the infinite Bethe lattice 
for every vertex and it is close to 1 for vertices belonging to a finite subset of 
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the Bethe lattice which are far away from the "boundary" (see later on for 
precise statements). In our case, we show that for an infinite random tree the 
characteristic ratio is a well-defined random variable, uniquely determined by 
a stochastic recursion. The transfer matrices will then contain elements with 
that distribution. We also consider deterministic trees that are strict subsets 
of the binary tree where the characteristic ratio can be computed explicitly. 
Next, we prove the exponential decay of correlation of height variables (as 
in the Bethe lattice case), and show that for a branching probability p suf- 
ficiently small, but still supercritical (p > 1/2), i.e., the branching process 
survives with positive probability, avalanche sizes decay exponentially. This 
shows a transition between exponential decay of avalanche sizes, for p small, 
and power law decay for p close to (possibly only equal to) one. 

Our paper is organized as follows. First we recall some basic material 
about the ASM on trees and the recursive technique developed in |5] . Second 
we study the recursion for the characteristic ratio and show it has a unique 
solution for the random binary tree. Finally we give quenched and annealed 
estimates of the eigenvalues of the product of n random matrices, which we 
apply in the study of correlation of height variables and avalanche sizes. 

2 Abelian sandpile model on subtrees of the 
full binary tree 

We summarize here the basic and standard objects of the abelian sandpile 
model on a (general) tree. More details can be found e.g. in [5l ITT] . 

2.1 Rooted and unrooted random trees 

We denote by the rooted binary tree of n generations, and by the 
rooted infinite binary tree. For a more general tree ^ we write if we 
want to indicate that the tree has root i. The rootless infinite binary tree or 
Bethe lattice is then obtained by joining two infinite rooted binary trees by 
a single edge connecting their roots. 

A random binary tree of generations with branching probability p G 
[0, 1] is a random subset ^ of obtained as follows. Starting from the 
root, we add two new vertices, each connected with a single edge to the root 
(resp. no vertices), with probability p (resp. (1 — p)), and we iterate this from 
every new vertex independently for A^ generations. By letting A^ — i- oo we 
obtain the full binary branching process. Joining two independent infinite 
copies of this process by a single edge connecting their roots creates the 
rootless random binary tree. This last procedure is of course identical to 
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create the non-random rootless binary tree from non-random rooted binary 
trees. 

2.2 Height configurations and legal topplings 

For ^ a finite subtree of the Bethe lattice, height configurations on ^ are 

elements r] G {1,2,.. .}'^ := ^^.y. For r\ e and u G r]u denotes the 
height at vertex u. A height configuration ri G J^?^ is stable if r]u G {Ij 2, 3} 
for all M G Stable configurations are collected in the set = {1, 2, 3}"^. 
For a configuration rj G we define the toppling operator via 

where A is the toppling matrix, indexed by vertices u,v & ^ and defined by 

A„„ = 3, Auv = —1 if u, V are neighbors in £^ (1) 

{u,v neighbors in ^ means that an edge of ^ connects u to v). In words, 
in a toppling at u, 3 grains are removed from u, and every neighbor of u 
receives one grain. 

A topphng at u G ^ is called legal if 77^ > 3. A sequence of legal topplings 
is a composition T^^ o . . . o (rj) such that for all /c = 1, . . . , n the toppling at 
Uk is legal in T„j._^ o . . . oTu^{ri). The stabilization of a configuration 77 G 
is defined as the unique stable configuration S^{rj) G that arises from 77 
by a sequence of legal topplings. 

2.3 Addition operator and Markovian dynamics 

For ^ a finite subtree of the Bethe lattice, and for m G the addition 
operator is the map '■ ^ sr ^ ^ s defined via 

a^r]^y{r] + 5u) (2) 

where 5„ G {0, 1}^ is such that 5u{u) — 1 and 5u{z) — for z e , z u. 
In other words, a^ri is the effect of an addition of a single grain at in 77, 
followed by stabilization. 

The addition operators commute, i.e., a^a^ = a^a„. This is the well- 
known and crucial abelian property of the sandpile model. 

The dynamics of the sandpile model is then the discrete-time Markov 
chain {77(71), ti G N} on defined via 
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(3) 



where Xi are i.i.d. uniformly chosen vertices of 

Given a stable height configuration rj and m G we define the avalanche 
Av{u, rj) induced by addition at u in 77 to be the set of vertices in ^ that 
have to be toppled in the course of the stabilization of r] + 6u- 

2.4 Recurrent configurations and stationary measure 

The recurrent configurations of the sandpile model form a subset of the 
stable configurations defined as follows. A configuration rj G Q,3r contains a 
forbidden subconfiguration (FSC) if there exists a subset S G ^ such that 
for all u G S, the height in u is less than or equal to the number of neighbors 
of u in S. The restriction of rj to S is then called a FSC. A configuration 
is allowed if and only if it does not contain a FSC. Recurrent configurations 
coincide with allowed ones, and are collected in the set ^^r. 

We denote by =^^(^2,57) the set of probability measures on The Markov 
chain ^ has a unique stationary probability measure /^^ G ^(fi^) which 
is the uniform measure on the set 



where 5^ is the point mass concentrated on the configuration r]. 

2.5 Specific properties of the sandpile model on a tree 

In this section, for the sake of self-containdness, we briefly summarize some 
basic facts from the paper [5] which we need later on. In a subset of the 
Bethe lattice, the distance between two vertices is defined as the length (i.e., 
the number of edges) of the shortest path joining them. A vertex is a surface 
vertex if it has a number of neighbors strictly less than 3. 

2.5.1 Weakly and strongly allowed subconfigurations 

The class of allowed configurations can be divided into weakly and strongly 
allowed ones. Let T be a rooted finite tree with root u and extend it with one 
vertex v and one edge < u,v >. Consider an allowed configuration ^ on T. 
We put height 1 at vertex v and investigate the derived configuration ^' on TU 
{v} such that ^ is the restriction of ^' to T, and = 1. If ^' has no FSC, then 
we call ^ strongly allowed on T, otherwise weakly allowed. We give in Figure 
1 an example of weakly and strongly allowed configurations. On one hand, 
if rju = 2, then there exists a forbidden subconfiguration on S* = {v,u,ii}. 
On the other hand, if rju = 3, there are no forbidden subconfigurations. 
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Figure 1: Example of a tree with root u 



2.5.2 Characteristic ratio and recursion 

Let <^ be a finite tree, rooted or not. A key quantity in the analysis of [5] 
is the characteristic ratio x{£^) G [1/2, 1] of ^ between weakly and strongly 
allowed configurations. For the empty tree ^ = we put x{0) = 0. 

For an infinite tree =^ we say that the characteristic ratio is well-defined if 
the limit Ymva^-i^^'^ x{^') exists, where \im.f^i^y is taken along the net of finite 
subtrees. That is, \im,yi^,^ x{^') = a means that for every e > 0, there exists 
<5q' C a finite subtree of £^ such that, for all finite subtrees D we 
have \x{£^') — a\ < e. 

The characteristic ratio satisfies a recursion property both for rooted and 
unrooted trees: 

^ ' 2 + x(Ti) + x(T2) 

where T^, are two non-intersecting subtrees of the tree T defined as follows. 
If T = T° is a rooted tree with root o, then T^^T"^ are obtained by deleting 
the root o and splitting the tree into two subtrees whose roots are the two 
descendants of the root o, see Figure 2. If T is an unrooted tree, we pick one 





Figure 2: Example of a tree with root o which splits into and 
of its edges, < m, f >, we delete this edge and split T into two rooted subtrees 



and whose roots are u and v, see Figure 3. This recursion (ll) holds for 



Figure 3: Example of a rootless tree which splits into and 

finite trees, and by passing to the limit for infinite trees, provided the limits 
defining x(T), x(T^), x(T^) exist. Dhar & Majumdar obtain from Q that 
when the distance from the root of T to the nearest surface vertex tends to 
infinity (what they call "deep in the lattice"), the characteristic ratio x{T) 
tends to 1. The authors then derive several explicit quantities (such as height 
probabilities of a vertex u deep in the lattice) by replacing the characteristic 
ratio by 1. 

2.5.3 Transfer matrix approach 

The transfer matrix approach allows to compute the two-point correlation 
functions and avalanche size distribution. Let u, v be two vertices in the tree 
at mutual distance n. They determine n + 3 subtrees Ti, . . . Tn+3, see Figure 
4. The number of allowed configurations when fixing the heights at u and v 
can then be obtained via a product ^ of n + 3 two by two matrices with 
elements determined by the characteristic ratios of the subtrees Ti, . . . , Tn+3 
(see Section 4 below for the precise form of these matrices). In particular 





Figure 4: Example of a tree with n + 3 subtrees 
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for the Bethe lattice, all the matrices involved in this product are equal to 
2 2 



, because the characteristic ratios of all the involved subtrees are 

-L o 

equal to one. Let ^ he a finite or infinite subtree of the full binary tree. The 
two-point correlation function, i.e., the probability that two vertices u,v at 
mutual distance n have height i resp. j, for i,j G {1,2,3}, is equal to (see 
[3 Section 5, eq. (5.11)]) 



where A_(^)/A_|_(^) is the ratio of the smallest and largest eigenvalues of 
the matrix ^ and aij are some numerical constants depending on For 
the binary tree, A_(^) = 1 and A+(^) = 4"^. 

The avalanche size distribution is determined by the inverse of X+{^). 
For the binary tree, upon addition of a grain at a vertex u, the probability 
that the avalanche Av{u,ri) is a given connected subset ^ of cardinality n 
containing u is equal to 

/i^(M«,r/)=^) = ^, (6) 

for some constant C, independent of the shape of the subset Since there 
are 4"n~^''^(l + o(l)) connected subsets of the Bethe lattice of cardinality n 
containing u, one concludes that for large n (see [5i Section 6, eq. (6.13), 
(6.14)]), 

^,y{\Av{u,r])\ =n) ^ rT^^'^ (7) 
i.e., the tail of the avalanche size distribution decays like rr'^l'^ . 

3 Some characteristic ratios 

On the Bethe lattice x(^*) = 1 for the (infinite) rooted subtrees (i = 1, 2, 3) 
attached to every vertex since every is an infinite rooted binary tree (see 



Subsection 2.5.2). This property considerably simplifies the analysis of |5] 



and is no longer valid in the inhomogeneous or random cases studied here. 

3.1 The characteristic ratio of the random binary tree 

Let be the random binary tree of n generations starting from a single 



individual (the "zero-th" generation) at time n = (cf. Section [2T). Then 



for n > 0, x{£^n) satisfies the recursive identity [5l Section 3, eq. (3.12)] 



8 



where 

■'^ ' ^ 2 + u + v \l + u 1 + vJ ^ ' 

and are the (possibly empty) subtrees emerging from the (possibly 

absent) individuals of the first generation. For a tree ^ consisting of a 
single point we have 

= I (10) 
because heights 2, 3 are strongly allowed and height 1 is weakly allowed (cf. 



Section 2.5.1). This value 1/2 can also be obtained from the recursion ([8]) 
by viewing a single point as connected to two empty subtrees (for which 
x{o) = 0). Notice that if u,v E [0, 1], then f{u,v) < 1. Therefore, we view / 
as a function from [0, 1]^ onto [1/2, 1]. 

Lemma 3.1. For every finite subtree 3^ of the Bethe lattice, the ratio x{^) G 
[1/2,1]. Moreover, on [0,1]^ the function f defined by ^ is symmetric, i.e., 
f{u,v) = f{v,u), and increasing in u and v, i.e., forui < U2, 

/(mi, v) < f{u2, v), for all v (11) 

and analogously in the other argument. 

Proof. The proof is straightforward and left to the reader. □ 



Proposition 3.1. There exists a random variable X^o such that x{^n) — ^ 
in distribution as n ^ 00. 

Proof. Denote by /i„ G the distribution of x(e5^). The recursion 

([s]) induces a corresponding recursion on the distributions /in+i = ^(/^n)- We 
show that ^ is a contraction on the set ^{[1/2, 1]) of probability measures 
on [1/2, 1] endowed with the Wasserstein distance. This implies that it has 
a unique fixed point /i* and from every initial /io, yUn — ^ A** in Wasserstein 
distance and thus weakly. 

Let g he a. Lipschitz function. Then we have 

J gix{X))f-in{d{x{X))) = 5^(?(x(r))/i„(x(5;) =x(T)) 

«e{o,i,2} 

+ 9{x{T))iin{x{:^n) = x{T)) (12) 
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T2 




Figure 5: The trees T^,T^,T^ 



where T^,T^,T'^ are the three trees which cannot be spht into two subtrees 
both non-empty, see Figure 5. Every other tree T (appearing in the last sum 



of (12)) can be spht into two subtrees both not reduced to a single point. 



Then using (\8h, the expression in (12) becomes 



9\ ^)(l-p)+^7o/(^,^ 



+ 2p\l-p) [gof(x{X-i 



1 



yu„_i((i(x(5;_i))) 



+ / [gof{x{X-i),x{X-i)) ]f^n~Mx{X-i))) 



Lemma 3.2. The function ^ on ^([1/2, 1]) defined by 
g{x)^{lj){dx) 

+ 2p^{l-p) j go f(^x,^fi{dx 



9\ l)(^-P)+9(l 




9 o fix,y)n{dx)n{dy) 



(13) 



(14) 



is a contraction on ^{[1/2,1]) endowed with the Wasserstein distance. The 
contraction factor is hounded from above by 8/9. 

Proof. Denote by ^ the set of Lipschitz functions 5^ : [1/2, 1] — )■ M with 
Lipschitz constant less than or equal to one, i.e., such that \g{x) — g{ii)\ < 
\x — y\ for all x,y. We use the following two formulas for the Wasserstein 
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distance of two elements /x, z/ of ^([1/2, 1]) [6]: 



d{fi, p) = sup 



gdfi — / gdu 



d{fi, v) =ini{ / \x- y\¥{dxdy) : = /x, P2 = 



(15) 



(16) 



where in the last right hand site the infimum is over all couplings P with first 
marginal Pi (resp. second marginal P2) equal to yU (resp. v). 



To estimate d{^ (fi) , ^ (u)) we start with the first formula (15) 
d(^(/i),^(i/)) = sup [ g{x)^{fi){dx)- [ g{y)^{u){dy) 



< 2p^(l — p) sup 



9 ° f{x,-)fJ.idx) - / g o fiy,-)u{dy) 



+ p sup 




g o /(x, x')^{dx)fi{dx') 




9° f{y,y')^{dy)u{dy') 



;i7) 



=: 2p^{l - p) sup A{ii, V, g) + p^ sup B{^, v, g). 
Now use the definition of / and the fact that x,y E [1/2, 1] to estimate 



f[y 



< -\x — y\. 

This gives, using the Lipschitz property of g and a coupling P of /i and z/: 



9° f[x,^ ]fi{dx) 



9° f[y,^ ]^{dy) 



< 



< 



< 



9° fix, 



1 



9° f{y, 



1 



F{dx, dy) 



1 



f[y 



X, dx' , dy, dy') 



4 



\x — y\F{dx, dy). 



Taking now the infimum over all couplings P, using (16) to bound ( |18[ ), and 
taking the supremum over g & ^ yields 



snpA{fx,u,g) < -d{n,u). 
ge^ 4 



(19) 
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To estimate the term B{n, u, g), use the elementary bound (since x, x', y, y' G 
[1/2,1]): ^ 

\f{x,x) - f{y,y')\ < -{4\x-y\ +A\x' -y'\). 

We then have, using the Lipschitz property of g, and a couphng P' = 
F'{dxdx'] dydy') of ® and v ® v 




g o /(x, x')ji{dx)jj,{dx') 




9° f{y,y')^{dy)u{dy') 




\f{x,x') - f{y,y')\F'{dxdx';dydy') 

< ^ J y\F'{dxdx'] dydy') + ^ J \x' - y'\F'{dxdx'; dydy'). 

Taking the supremum over g ^ and infimum over the couphngs P', we 
find 



sup B{n,iy,g) < -d{^i,v). 



(20) 



Combining (19), (20) with (17) we arrive at 



d{^{^^),^{v)) < [^-p\l-p) + -p^j d{^^,v) < -d{f^,u) 
where in the final inequality we used the elementary bound 



-/(I -p) + -p^ =p 



1 7p\ 8 



2 ' 18 / - 9 



□ 



This proof of Lemma 3.2 completes the proof of Proposition 3.1 □ 



Remark 3.1. The random binomial tree is such that every vertex has two 
children with probability p^ , 1 child with probability 2p{l — p) and children 
with probability (1 —p)"^, p G [0, 1]. For this case the same idea as in Lemma 



3.2 gives a recursion leading to a contraction in the Wasserstein distance and 
hence x{,%i) has also a unique limiting distribution. 



3.2 The characteristic ratio of some deterministic trees 



The recursion (|8j) allows also to compute the characteristic ratio for certain 
(deterministic) infinite subsets of the full binary tree in terms of iterations 
of sections of /. 
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3.2.1 Infinite branch 

First, consider a single branch of length n, i.e., the tree consisting of 

a root and n > 1 generations of two individuals each (see Figure 6). Using 




Figure 6: Example of a tree with a single branch of n generations 



(10) we obtain the recursion 



^ (^nT") =f[\^ ^r"" ) (21) 



which by Lemma 3.1 gives that x(<^™") is monotonically increasing in n 
with limit 

x^=x ) = (22) 

which is the unique positive solution of 

x = /Q,x^. (23) 

Remark 3.2. We can generalize this "backbone" tree to a general "backbone- 
like" tree where at each point the same finite tree T is attached (it is a single- 
ton in the backbone case). The characteristic ratio equals the positive solution 
of the fixed-point equation 

X* = f{x{T),x*) (24) 
where x{T) denotes the characteristic ratio of the tree T, i.e., 

X* = ^(-1 + V5 + 4x(T)). (25) 
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3.2.2 Finite perturbations of a single branch 



Attaching a finite subtree T of the full binary tree at level n in the 
infinite single branch =5^*" leads to a tree ^^'^'^ with characteristic ratio 



(26) 



where (p{x) = /(l/2,x) (cf. (23)) is applied n times. This shows that 
inserting a finite tree T at level n has an effect on the characteristic ratio 
that vanishes in the limit n — )■ oo, exponentially fast in n. 

Moreover, since x{T) > 1/2 and x t-?- f{x,y) is monotone for all y (by 
Lemma 3.1), we have from (26) 



x{^r)>^o<p---o^{f{7;,<^::n 



x(5; 



min\ 
oo ) 



(27) 



From ( p7| ) and (22) we conclude that for every infinite subtree =^ C for 
which x{£^oo) exists, 

= Xi^^n < Xi^oo) < X(^oo) = 1. 



4 Transfer matrix and eigenvalues: uniform 
estimates 

In the analysis of the two point correlation function and of the avalanche size 
distribution, one is confronted with the problem of estimating the minimal 
and maximal eigenvalues (denoted by A_(^) and A_|_(^)) of a product of 
matrices of the form 

n 

^(xi,...,x„) = JJm(x,) (28) 

with 



=1 



M(..)=(l + -' ^ + ^;) (29) 

where the the characteristic ratios of some recursively defined sub- 

trees. 

When Xi = 1 for all i, this exactly corresponds to the analysis in [5] 



see the above Subsection 2.5.3). More precisely, in that case, one needs the 



minimal and maximal eigenvalues (denoted by A_ and A+) of 

^(1,1...,1)= ( 1 3) 
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which are A_ = 1 and A+ = 4". This leads to a decay of the covariance 
proportional to 1/4" = A_/A+ and decay of avalanche size asymptotically 
proportional to 

4n J ^3/2 

where 

Ar,= \{'€ d SB^:'^ connected, |^| = n, ^ 9 o}| (30) 

denotes the cardinality of the set of connected clusters of size n containing the 
origin o. In the general case the decay of covariance can be estimated by the 
ratio A_(^)/A+(^), and for the decay of the avalanche size distribution, 
one needs to estimate A+(^) as well as the analogue of A„. In the case 
of a branching process, the Xj's appearing in the matrix ^(xi, . . . ,a;„) are 
independent random variables with distribution /i* defined in the proof of 



Proposition 3.1 



In this section we therefore concentrate on the estimation of the eigen- 
values of a matrix of the form ^( ) for general Xj's. 

Lemma 4.1. 1. For all n and all Xi [1/2,1], (1 < i < n), the eigenval- 
ues of ^{xi, . . . , Xn) are non-negative. 



2. We have the inequality 



3. We have 

A_(^) det(^(a;i, . . . , Xn)) 



^-^^^<c(i] . (31) 



(1 + 0(1)) (32) 



Tr{J({xx, ■ ■ . ,x„))^ 
where o(l) tends to zero as n ^ oo, uniformly in the choice of the Xi 's. 

Proof. The eigenvalues are given by 



^{a± Va^ - 46) (33) 



' 2 

with a = Tr(^) and b = det(^) (we abbreviate ^ for ^{xi, . . . ,x„) in 
these expressions). For n = l, a^ — 46 = 5 + 4x>0. Hence 

3 + 2x ± V5 + 4x 
A± = ^ • (34) 

For n > 2, we estimate, as in [TT] 

n 

a = Tr(^) > JJ(2 + a;i) 

i=l 
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and using 



we have 



n 



1=1 ^ 



i=l 

-2 



< 



(35) 



(36) 



In particular, 46 < which imphes that the eigenvalues are real and non- 
negative (a > 0). Inequality (31) then follows from (36), (33). 

Given that (by (36)) h = h/a? tends to zero as n — )■ oo at a speed at least 
C(4/9)", we have 



A_(^) _ 1 - VI -4/i 
A+(^) ~ 1 + Vl + 4/1 

which proves the third statement. □ 



h{l + 0{h)) 



Lemma 4.1 shows that the ratio A_(^)/A+(^) behaves in leading order 
as det(^)/(Tr(^)^). To estimate this ratio, we need to estimate Tr(^) 
from below. We start with a useful representation of ^(xi, . . . , Xn)- 



Lemma 4.2. Define 



1 1 
1 





1 1 



a=(cn,...,a„)e{0,l}" i=l 



Then we have 

Jl[Xx, . . . ,Xn) 

where yi = {1 + Xi) and where 

n 

= n [ETEl--^] . 

i=l 

Proof. We have 

n 

M (Xi, . . . , X„) = JJ(l/i^l + E2). 



a] 



(37) 



(38) 



(39) 



2 = 1 



The result then follows from expansion of this product. □ 
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Lemma 4.3. 1. For all n>l, 



1 n 
1 



rpn 771 rpn rp 

nj2 — -C'2, -C'l -£^2 



n n 
1 1 



(40) 



(41) 



2. For all r >1, ki, . . . , kr, kr+i > 0, 

Tr(f[{E'l^E,)]=l[il + h) 

and 

Tr I (fliEl^E^) ] e1^-^^ ) = I 11(1 + h) ]{l + k, + (42) 



,i=l 



1=1 



,i=l 



A=2 



Proof. Identity (42) follows from (41) and invariance of the trace under 



cyclic permutations. To prove (41), use the expression in (40) for E1E2, and 



estimate the diagonal elements of the product 



ki h \ [ k2 k2 
1 1 J I 1 1 



which implies the result. □ 



1 1 



.1=2 



ki ki 
1 1 



Proposition 4.1. For all xi, . . . ,a;„ G [1/2, 1] we have 
Tr{^) > (f[{l + 2yi)] 2-"(i^/25)_ 



(43) 



.1=1 



As a consequence we have the following uniform upper bound 

A_(^) 



where 



7 



A+(^) 
4 



< C7" 



(44) 



25 



232/25 ^ 0.38854 



Proof. Remember definition (39) of <#"(«) for a = (ai, . . . , a„) G {0, l}". 



Using Lemma 4.3 we see that 



Tr(cf (a)) > 2"^(") 
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where ,J\^{a) = Y17=i '^ii^ ~ Q^i+i); - with the convention On+i = 1 - is the 
number of intervals of successive I's in the configuration a. Hence by (38) 
we obtain the lower bound 

Tr(^(a;i,...,x„)) = ) Tr(^(a)) 

a=(ai,...,a„)G{0,l}" \i=l / 

- ( nz/r ) 2^"=i"'(^-"-+i). (45) 

a=(ai,...,a„)e{0,l}" \i=l / 

Next since Ea=(ai,...,a„)G{o,i}" (nr=i(2y*)"0 = nr=i(l + 2l/i) =-Z,we rewrite 
ill ) 2^^=' = ZE, (2' (46) 

a=(ai,...,a„)e{0,l}" \i=l / 

where 

iEa(^(«)) = ^ ( n(2?/^)" ) ^(«) 

a=(«i,...,a„)e{0,l}" \i=l / 

defines a (|/j)i dependent probability measure on the a's. Now apply Jensen's 
inequality in (46) to obtain 

a=(oi,...,a„)e{0,l}" V«=l / 



Finally, use j/j = 1 + Xj G [3/2, 2] to estimate 

n 

^ (l + 2|/,)(l + 2y,+i) ^ 25 
which gives the following uniform lower bound for Tr(^): 

^Tr(^(xi,...,xO)>2-"(i6/2^). 
Combining this with det(^) = YYi=i ^1 obtain 

A-(^) <(i + o(l))2"(32/25)-Q^_ 2/^ 
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Finally we use Ui = (1 + Xj) with Xj G 
x)(3 + 2x)~^ is increasing to estimate 



[1/2, 1] and the fact that x i-)- (1 + 



2y^ 



< 



25 



which implies (44). □ 



5 Transfer matrix: annealed estimates 

In this section we look at the eigenvalues of ^{xi, . . . where now the 
Xj's are i.i.d. with a law /i on [1/2, 1]. We denote by P the joint law of the 
Xj's and by E the corresponding expectation. 
We start with the following lemma. 

Lemma 5.1. For all '-f > the eigenvalues of 



Ml) 

given by 



1 1 
1 1 + 7 



A,„) . ^^ + l^„v^ (47) 



are non-negative. 

Proof. Elementary computation. □ 



Theorem 5.1. Let A-i-(^) denote the largest, resp. smallest, eigenvalues 
of ^{xi, . . . ,Xn) where the Xi's are i.i.d. with a law supported on [1/2,1]. 
Denote 

r„.lU'^). ,48) 



Then we have: 

1. Concentration property: there exists C > such that for all e > 

P((F„-E(F„))>£)<e-^^'". (49) 

2. The limits 

L± = lim -logA±(^) (50) 

n— >oo n 

exist and satisfy L± = E,{L±) almost surely; moreover 

+ = 2E(log(l + Xi)). (51) 
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3. Upper hound 



lim r„ = L+ - L_ = lim E(F„) < log 



A- (7) 



(52) 



where A± are given by (47) and 'j = (E((l + xi) ^)) E [3/2,2]. 



Proof. We have, using (38) 



y = ^°g( E (nr-MTr(^(«)) I =:^(z/i,---,2/n). (53) 

,a6{0,l}" \i=l / 



Using that the weights 



wia] 



nr=ii/r"^Tr(^(«)) 



E.e{o,i}"(nr=iz/r'')Tr(^(«)) 



are non-negative and sum up to one, we compute the variation 



Li^p = sup iV-d/i,. , 



yi,...,bi,Vn 



< 



\ ae{0,l 



< — sup log 

n a,6e[l/2,l] 



w[a\ 



Vi 

log 2 

n 



.Vn) 



(54) 



Statement 1 is then an application of the Azuma-Hoeffding inequality. The 
first part of Statement 2 follows from Oseledec's ergodic Theorem [13], to- 



gether with (35) and the law of large numbers. To prove Statement 3, start 



from (38), (32), and use Jensen's inequality, combined with the mutual in- 
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dependence of = 1 + Xj, 1 < z < n, to obtain 



lim E( — lo; 



A+(^) 



2n ^\\_{J^) 



< 



lim El -log J2 (nyrMTr(^(")) 

y ae{o,i}" \j=i / 

\ae|0,l|" \j=i / 



\a6{0,l}" \i=l 



Tr(^(a)) 



log 



A- (7) 



(55) 



with = E(yf and A±(7) given by fi/H). Because of p5|), we finally 
derive (51) by the law of large numbers. □ 



6 Covariance and avalanche sizes 



6.1 Quenched and annealed covariance 

Theorem 6.1. Let ^ be a finite or infinite subtree of the full binary tree SS^. 
As before, denote by [i^ the uniform measure on the recurrent configurations 
Mg- of the sandpile model on ^ and let u,v & ^ be at mutual distance n. 
Then we have the following estimate for the two point correlation function 



\Cov{u,v,,^)\ :-- 



<C7". (56) 



In particular, Cov{u,v, ^) is absolutely summable, uniformly in the choice 
of the subtree 

sup ^\Cov{u,v,£^)\ < 00. (57) 



Proof. The first statement (56) follows from the expression of the covari- 
ance in Subsection 2.5.3, eq. ^5|, together with Proposition 4.1 Since the 
number of vertices at distance n in ^ from the origin is bounded from above 



by 2", and 7 < 1/2, (57) follows from (56). □ 
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Theorem 6.2. Let 3^[o,n] be the stationary binary branching process with 
reproduction probability p, conditioned to have a path from its root o to a 
vertex at distance n. Let fi,sr[o,n] denote the uniform measure on recurrent 
configurations on ^[o,n]. Let Cov{o,v{n), ^[o,n\) denote the covariance of 
the height variables at a and at a vertex v{n) at distance n from the root a. 
Then we have the following annealed lower bound on the covariance. 



lim l(-logE(|Coi;(o,t;H,^[o,n]))|) <logf^44 



< 



4-^7 



(58) 

where 7 ^ = E(?/]^ ^) and yi = 1 + xi, with xi distributed according to the 



measure of Proposition 3.1, and A± given in (47) 



Proof. The result follows from the expression of the covariance in Subsec- 
tion [2^ eq. (|5|), together with (pi). □ 



6.2 Avalanche sizes 

We start by defining a matrix associated to an avalanche cluster. Roughly 
speaking, the probability that the avalanche coincides with the cluster is the 
inverse of the maximal eigenvalue of this matrix. Let ^ C be a finite or 
infinite subtree of the full rootless binary tree, containing the origin 0. Let 
^ be a connected subset of ^ containing the origin. 

Definition 6.1. Let \^\ = n. The matrix associated to ^ in S/' , denoted 
by M(^) is defined as follows. To the vertices in are associated n + 2 
subtrees Ti, . . . ,T„+2 with corresponding characteristic ratios x(Tj). Then 

M,^,^n(l + f-) l + ^JI). (59) 

i=l ^ ^ 

Lemma 6.1. As before, for a stable height configuration rj, let Av{o, r]) denote 
the avalanche caused by addition of a single grain at the origin in rj, and fi^r 
the uniform measure on recurrent configurations on !J . Then there exist 
constants €±,€2 > such that 

ci {X4M{^)))-' < fi.^{Av{o, r/) = ^) < C2(A+ {M{^)))-' (60) 



where M{^) is the matrix of (59) associated to ^ of Definition 6.1 



Proof. It follows from the expression for fi^y{Av{o,r]) = ^) in Subsection 
2X3| eq. (g. □ 
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Definition 6.2. We denote by A„(o, ^) the number of connected subsets of 
edges of ^ containing the origin o and of cardinality n . 

1. The growth rate is defined as 

fi;(o, ^) = limsup -logAn(o, 5^). (61) 

n— ^oo ri 

2. We define the averaged growth rate as 

R = limsup — logE(y4„(o, ^)). 

The growth rate gives the dominant exponential factor in the number 
An{o, E.g. if <^ is the full binary tree, k = log 4, since 

A^{o, = C4"n-3/2(l + 0(1)) (62) 



see also [S], Subsection 2.5.2, For a stationary binary branching process we 
have the upper bound 

E(A„(o, ^)) < Cp^'Anio, ^oo). (63) 



For an exact expression E(y4„(o, ^)) from which inequality (63) follows im 



mediately, we refer to the Appendix, Proposition 7.2 

Theorem 6.3. 1. There exists C > such that for all n > 1 we have 

fi^{\Av{o, r])\=n)< CAr,{o, ^)4-"2"(i6/25)_ ^54^ 

In particular if the tree has growth rate k,{<^) < ||log2, then the 
avalanche size decays exponentially. 

2. For the stationary binary branching tree with branching probability p, 
we have the estimate 

E(/i,^(|At;(o,r/)| = n)) < CT^^^m {^-^^^ . (65) 
In particular, for 

p < 0.54511... 

the averaged (over the realization of the tree) probability of an avalanche 
size n decays exponentially in n. 
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3. For the binomial branching tree we have 

¥.{^i,^{\Av{o,ri)\ = n)) < C2-"(i'^/2V, 

In particular, for 



(66) 



p < 2-16/25 ^ 0.641713... 

the averaged ( over the realization of the tree ) probability of an avalanche size 
n decays exponentially in n. 



Proof. The first result, i.e., (64), follows from Lemma 6.1 , and the estimate 
(43) which gives the following lower bound on the largest eigenvalue 

A+(M(^)) = Tr(M(^)(l + oil)) > C74"2-"(i6/25)_ ^g^^ 



The second result follows from Lemma 6.1, (67), (63), (62) and Proposition 
17.11 below. The solution of 



is given by 



P+VP ^ 2-16/25 

2 



p ^ 0.54511 



(6J 



(69) 



The l.h.s. of (68) as a function of p is monotone, which yields that for 
p < 0.54511 the expected number of avalanches of size n decays exponentially. 

The third statement follows from Lemma 6.1, (67), (63), (62) and Remark 
[O □ 



Remark 6.1. From Theorem we conclude that avalanche sizes decay 
exponentially forp small enough. We know that on the full binary tree (cor- 
responding to p = 1) avalanche sizes have a power law decay. It is a natural 
question whether the transition between exponential and power law decay oc- 
curs at some unique non-trivial value Pc G (1/2, 1). This question is related 
to the behavior of a a random walk on the full binary tree, killed upon exit- 
ing a random subtree. Would this random walk have a survival time with a 
finite exponential moment, then we would have Pc = 1. Notice that for the 
corresponding problem on the lattice 1^^ , this random walk does not have a 
finite exponential moment of its survival time, because the tail of the survival 
time is stretched exponential (Donsker-Varadhan tail g-c*'*'^*''^^' However, 
on the tree in the annealed case corresponding to Theorem 6^, the decay of 
the tail of the survival time is not straightforward. Formally taking the limit 
d oo in the Donsker-Varadhan tail suggests that on the tree, the survival 
time of this random walk has a finite exponential moment, which points into 
the direction pc = I. See also for trapped random walk on a tree. 
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7 Appendix: The expected number of clus- 
ters containing the origin in a random bi- 
nary tree 

Proposition 7.1. Let p be the probability that a vertex has 2 children. Fur- 
thermore let An{o, denote the number of connected clusters of n edges 
containing o, given a realization of a tree with root o. Then there exists 
a constant C > such that for n > 1 we have 

E{A,,{o,^°))<C{n + l)r(^^^) . (70) 



Remark 7.1. For the binomial branching tree, where every vertex has 2 
children with probability p'^ and 1 child with probability 2p{l—p) (see Remark 



3.1), we have the upper bound 

E{An{o, sr°)) < crp"" (71) 

which follows from the observation that the number of n vertex animals con- 
taining the root of the full binary tree is bounded from above by 4", and in the 
case of binary branching tree, each of these vertices is present with probability 
p independently, which gives the factor p"' . 

Proof. If we define a„ to be the expected number of connected clusters 
containing o and with n vertices, then a„ satisfies the recursion relation 

n-l 

an = ]l{n=l} + P]l{n>2} ^ aiUn-i-l, U > 1. (72) 

1=0 

Furthermore, by definition we put ag = 1. Thus, going from vertices to edges, 

E(A„(o, ^°)) = an+i. (73) 
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Introduce then the associated generating function: 

oo oo n— 1 



n=0 



n=2 1=0 
\2/ 



(74) 



For X G 



0, 



-p + Vp 



[ '2p(l-p)J 
it is given by 

A{x) -- 



= 1 + X + px{A\x) - 1). 
, this power series is convergent and, using A{0) = 1, 



— ( 1 - a/1 - Apx(l + x(l - p))]. 
px \ J 



for a; G ( 



-p + 

M^-p). 



2p 



. Use (for z such that < 1) 



(75) 



^-^ \ n — 1 In 

n=l ^ / 



(76) 



to obtain 



oo n+1 



•4W = EE 

n=0 j=0 



2n\ fn + 1 



n 



J 



[i-py 



p" 



n + 1 



-X 



n+j 



(77) 



We put k = n + j, then n = k — j and j < n + 1 yields j < [{k + 1)/2J , 
hence the expected number of clusters containing the origin is equal to 



where 



L(fc+1)/2J 
j=0 



2{k - j)\ fk - j + 1\ f 1 - p 



(7J 



k-j 



P 



k-j + 1 



[k~m{k-2j + i)\ 



(9 



(79) 



with c = (1 — j9)/j9. All terms in the sum ( [78| ) are exponentially large in /c, 
hence the exponential growth of the sum is determined by its maximal term. 
Using Stirling's approximation n\ n^e~^^2'nn for the right hand side of 



(79), we obtain the upper bound 

e 



(k-j) 



7rV2 3'{k-23f-^o 
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^0) 



Define, for x e [0, (A; + l)/2], 



-X (k X 



1 k—x 



TT 



^/2 a;^(fc-2x) 



k-2x ' 



Tflis function ip attains its maximum at x = k{l — ^)/2, wliich combined 
witli (80), implies 



V2 \ 2^ 



(81) 



Plugging (81) into (78) and bounding from above induces (70). □ 



In the following proposition we give an exact formula for E(y4„(o, ^°)). 
Proposition 7.2. We have the following identity, 



E(A„(o, ^°)) 

/2(n + l^ 



/n + 2 n + 1 1 -{1-P )\ (82) 

n + 2\ n + 1 / \ 2 2 2 p 

where 2-Fi(", ■) denotes the hypergeometric function defined as 

(a)„(6)„2;" 



2Fi{a,b,c,z) = ^ 



n=0 



c n n! 



^3) 



and where the Pochhammer symbol (a)„ is defined by 



ifn = 0, 
a{a + l)...(a + n — 1) if n > 0. 



M) 



Proof. Let us first remark that for \z\ < 1, real a,b and c 7^ — m where 



m G N, the series in (83) is well defined. In our case z = —(1 — p)/p, 



a = —{n + l)/2, b = 1/2 — (r2 + 1) and c = (« + 1) — 1/2 ( see also [IJ for more 



details about the hypergeometric function). The claim (82) follows from (73) 



and (78) once we show the identity 



2(A;-j)\ A-j + 1 



L(fc+1)/2J 

E 

j=0 

1 f2k\ f k + 1 k 1 
~ITi[k r \ 2~'~2'2 



1 — p 
p 

-k'- 



k-j + 1 

i-pY 



^5) 



p 
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If a or 6 are negative integers, using (84) we see that the series (83) defining 
2-Fi is a finite sum. Therefore, 



k+1 k 1 , 1-p 

k, 

P 



2 ' 2' 2 



j=0 



We have 



r(a + n) 



r(a) • 

Next we use the functional identities for the Gamma function (see [T6] ) 

— 7r 



and recurrence relation 



to rewrite 



z sin(7r2;) 
zT{z) = T{z + l) 

(-¥), (-1), 



in terms of Gamma functions with positive arguments. Thus 

k + i\ _ r(-^+j) 
r(-¥) 

r + 1) sin (vr^) 



(86) 
(87) 

(88) 
(89) 
(90) 



(91) 



which gives 



with 



A 



and 



B 



r(¥ + i-^0 sin 

r(^ + i)r(| + i)r(fc-j + i) 
r(¥ + i-j)r(| + i-j)r(^ + ^) 

sin (tt (^)) sin (f ) sin (vr (fc - | - j)) 



(92) 



(93) 



(94) 



sin (vr - j)) sin (tt (| - j)) sin (vr (A; - i)) ' 

Furthermore expanding numerator and denominator of this expression gives 
B = 



-ly . To rewrite A, we use (89j) and the duplication formula 

r 



T{z)T {z + 



)l-2z 



7rr(2^) 



(95) 
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once for z = {k + l)/2: 



k + l . + 1 , 1\ ^k + l_. 



H^^VH^+2j-^^'^^r(* + « (96) 

and another time for 2; = (fc + l)/2 — j: 

r - ^' + 1) r - + 1) = v^(^ - ^''-"''nk - 2j + 1), 



0f2-^+2J-ir(A;-2j +2). 



(97) 



H6I1C6 

^"r(A; + l)^' ^ r(fc-j + 2) ^r(fc-2j + 2) ^^^^ 
and since the Catalan numbers satisfy 

1 /2k\ 4^ T(k + l) 



k + l\k) 0Fr(fc + 2) 
we can write 



v^r(A; + 2) j! 



A^-i T{k-i + \) T{k - J + 2) 

X 



(100) 



^ V{k-3 + 2) r(A; - 2j + 2)r(j + 1) 
which is equal to 

and yields the claim. □ 

Remark 7.2. For p = 1, we recover the classical formula for the number 
of animals of n edges containing the origin in a binary tree (cf. Subsection 



n + 2\n + l / n large 

for some constant C . 

Proof. This follows from the fact that for p = 1, 2Fi{a, b, c, 0) = 1. □ 
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